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Abstract 


A three-dimensional (3D), two-phase, isothermal model of direct methanol fuel cells (DMFCs) was employed to investigate effects of electron 
transport through the backing layer and the land in bipolar plates. It was found that the electronic resistance of the backing layer, affected by 
backing layer electronic conductivity, backing layer thickness and flow channel width, played a relatively important role in determining the current 
density distribution and cell performance. In order to ignore the electron transport effect on the average current density, the minimum electronic 
conductivity of the backing layer has to be 1000 S m~t, with the relative error in the average current density less than 5%, under the given conditions. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


AS a promising candidate for next-generation portable power 
sources, the direct methanol fuel cell (DMFC) has attracted 
tremendous research and industrial interests. However, in order 
to compete with lithium-ion batteries, a portable DMFC system 
must overcome several key technical challenges, i.e., low rate of 
methanol oxidation kinetics, methanol and water crossover, ther- 
mal management, and catalyst poisoning and loading issues. As 
aresult, there is an urgent need for understanding, prediction, and 
optimization of various interactive transport and electrochemical 
processes that occur in portable DMFCs. 

Focusing on either one or two dimensions, early DMFC mod- 
els [1—6] were developed to study the mass transport phenomena, 
electrochemical processes, and their interactions, without con- 
sidering two-phase effects, which is found to be of paramount 
importance to understand DMFC behaviors. Meyers and New- 
man [7—9] developed a theoretical framework that describes the 
equilibrium of multicomponent species in the membrane in a 
three-part series. However, the effect of flow and stoichiometry 
was not considered. 
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In recently developed 3D polymer electrolyte fuel cell 
(PEFC) models, computational fluid dynamics (CFD) methodol- 
ogy has been employed to integrate electrochemical processes 
with water/proton co-transport to enable multiphysics model- 
ing and large-scale simulation. However, compared with H2/air 
PEFCs, the 3D mathematical modeling has not been widely 
applied to DMFCs, due to more complex interactions between 
two-phase species transport and electrochemical reactions. Most 
recently, Liu and Wang [10] presented a 3D, two-phase DMFC 
model, based on multi-phase mixture (M?) modeling framework 
of Wang et al. [11—13], to reveal an intricate interaction between 
species transport and electrical current in three dimensions. This 
model has been extended to further study the net water trans- 
port coefficient distribution and interfacial liquid water coverage 
effect by Liu and Wang [14]. A more complete review of DMFC 
modeling in the literature is recently given by Wang [15]. 

Among the 3D DMFC models published to date, electron 
transport in the catalyst and backing layers, and in the bipolar 
plate is ignored by assuming a sufficiently large electronic con- 
ductivity and consequently constant electronic phase potential 
in these materials. This assumption is valid likely in the lands in 
the through-plane direction, of which the electronic conductiv- 
ity is of the order of 20,000 S m7! for graphite plates, but could 
be inappropriate in the backing and catalyst layers, in which 
the effective electronic conductivities range only from 300 to 
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Nomenclature 


ajo total exchange current density (A m7?) 
c molar concentration (mol m7?) 

D diffusivity (m? s7!) 

F Faraday’s constant (96485 C mol`!) 

i operating current density (A m~?) 

j volumetric current density (A m7?) 

K, 


€ reaction constant 
m mass source term in governing equations 
ndg electro-osmotic drag coefficient 
p pressure (Pa) 
Ri ohmic resistance (i= land, backing or channel) 
(Qm?) 
s liquid water saturation 
sk molar source term of k species 
u velocity (m s7!) 


U? standard equilibrium potential (V) 
Veell cell voltage (V) 
Z width (in-plane direction) (m) 


Greek letters 


Qi current transfer coefficient at anode or cathode 
(i=aorc) 

ô thickness (m) 

E porosity of porous medium 

@ phase potential (V) 

y advection correction factor 

Ni overpotential of anode or cathode (i =a or c) (V) 

K ionic conductivity of membrane (S m7!) 

p density (kg m~?) 

o electronic conductivity of electrode (S m7!) 

Subscripts 

a anode 

c cathode 

cata catalyst layer 

e electrolyte 

eff effective 

g gas phase 

int interface 

l liquid phase 

m membrane 

s solid or electronic phase 

xover crossover 

Superscripts 

MeOH methanol 

O2 oxygen 

k species 

ref reference 


500S m7! in the through-plane direction. In recent studies of 
modeling on PEFCs, the electron transport in PEFCs has been 
addressed [16-21]. It was concluded that by solving the electron 
transport equation in the catalyst layers, the gas diffusion lay- 


ers (GDL), and the current collecting plates on both anode and 
cathode sides, it enabled further capabilities in modeling fuel 
cells, including direct incorporation of the contact resistances 
between the membrane electrode assembly (MEA) and GDL, 
and between GDL and the current collecting land. 

In this paper, a 3D isothermal, electrochemical and transport 
fully coupled DMFC model is described, with focus on elec- 
tron transport. Additional to the complete set of conservation 
equations of mass, momentum, and species presented in [10], 
an electron transport equation is included in the catalyst and 
backing layers, and the land in bipolar plates in order to inves- 
tigate the effects of electronic resistance in the DMFC. At the 
same time, a proton transport equation is solved inside MEA to 
accurately account for electrolyte ionic resistance. The effects 
of the lateral electronic resistance on the current distribution and 
cell performance are studied in detail. 


2. Model development 


The two-phase steady-state DMFC model in this work con- 
sists of nonlinear, coupled partial differential equations of 
conservation of mass, momentum and species. These can be 
presented in concise form as follows [10]: 


dep) 


Continuity : J + V- (pu) = (1) 
1 fo 1 
Momentum: — | (ou) + -V - (puu) 
e| at E 
=-Vp+V-Tt+ Su] (2) 


i ack . 
Species : a +V. {yuc"} =V. [Di eVect + D$ eV Cel 


auig® (3) 


where p, u, and c* denote the two-phase mixture density, velocity 
vector and molar concentration of species k, respectively. The 
species considered here are methanol, water, and oxygen. The 
various two-phase properties, source terms, and thermophysical 
properties identified for various regions of a DMFC, as well as 
necessary boundary conditions, have been detailed in [10-13] 
and thus are not repeated here: 


Proton transport: O=V- (keff V Be) + So (4) 


j anode catalyst layer 
whereSo = 


— je + jxover cathode catalyst layer 


This equation indicates that protons are generated in the 
anode catalyst layer and consumed in the cathode catalyst layer. 
In the cathode catalyst layer source term, the methanol crossover 
flux through the membrane, jxover, is caused by electro-osmotic 
drag and diffusion, given by 


Jxover — y ( ,McOH i j DMeOH(eMeOR iat /m) (5) 
Gp O F Scata 
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In additional to proton transport, an electron transport equa- 
tion is employed in the current 3D model, which will be solved 
numerically in the backing and catalyst layers, as well as in cur- 
rent collecting lands on both anode and cathode sides in order 
to fully investigate its effect on the current distribution and cell 
performance: 


Electron transport: O= V - (Os eff V Bs) + So (6) 


where 


So = [ , i l 
Je — Jxover 

The source terms in the above equation correspond to the 
methanol oxidation and oxygen reduction reactions, creating or 
consuming electrons in either anode or cathode catalyst layer. 
Note also that electrons always flow from low potential to high 
potential location. 

Electrochemical kinetics in the proton and electron equations 
is expressed by Tafel approximations of Butler—Volmer equation 
in the anode and cathode, respectively, such that 


anode catalyst layer 
cathode catalyst layer 


ajrel MeOH |. ata expC (œa F/ RT Na) 


(7) 
MeOH | ata + Ke €xp((@aF/ RT )na) 


J= 


where na = Ds — e — U? with ®, =0 at anode current collect- 


ing land boundary. 
And 
02 
. . 2 f Cy? cata Ac F 
Je = J + Jxover = ajf ( Quel (1 — s)exp (-Fn.] 


(8) 


where no = Ø; — De — U? with ©, = Veen at cathode current 
collecting land boundary. 

As can be seen, besides intricate couplings between species 
transport equations, a numerical problem unique to 3D modeling 
of DMFC is the strong coupling between the species transport 
equations and electrochemical kinetics that are dependent on 
species concentrations. Therefore, the model must seek a stable 
and efficient solution procedure. 

In this paper, the present DMFC model is implemented into 
commercial software, Fluent, through user-defined functions, 
or UDFs, provided by Fluent in DEFINE macros, which are 
employed to implement sources terms, diffusion coefficients, 
material properties, and boundary conditions. The species equa- 
tions of methanol, water and oxygen, as well as the proton 
and electron transport equations, are implemented using user- 
defined scalars (UDS). More detailed numerical procedures have 
been described in [10]. 


3. Results and discussion 
3.1. Theoretical analysis 


Fig. 1 displays a schematic of the current collecting land and 
backing layer. The typical thickness of the land (djanq) is 2 mm, 


Stand Spacking 


Z 


‘channel 


Fig. 1. Schematic of the current collecting land and backing layer. 


thickness of the backing layer (Spacking) is 0.3 mm and the half- 
width of the flow channel (Zehanne) is 0.5 mm. If the effective 
electronic conductivities of the land in the bipolar plate (i.e. 
graphite) and backing layer (i.e. carbon paper) are assumed to be 
20,000 and 300 S m~!, respectively, the typical electronic resis- 
tance in the current collecting land, Riana, can thus be calculated 
as 


Oland 2mm 
Riana = = 


= = == 1 mQ cm? (9) 
Olana 20000 S m 


This indicates that the electron resistance in the land is very 
small, and therefore, always negligible. 

The ohmic resistance across the backing layer in the through- 
plane direction, Rpacking, is 


Sbacking = 0.3 mm 


Obacking 300 S m7! 


Rbacking = = 10m2 cm? (10) 
And, the electronic resistance of the backing layer in the in- 
plane direction, Rehannel, 1S 


Z channel E 0.5 mm 


= z = 16.7 M2 cm? (11) 
Obacking 300 S m7 


Rchannel = 


The electron phase potential variation is proportional to both 
electron resistance and average current density. Since the aver- 
age current densities in the through-plane and in-plane directions 
are different, the electronic resistances along the two directions, 
Rbacking and Rchannel, are not directly comparable. Assuming that 
the average current density under the gas channel is a constant i, 
and the average current density through the backing layer thick- 
ness is uniform, the lateral electronic phase potential variation 
in the backing layer can then be derived as 


Zchannel 1 Z. 
Z = channel 
A®, = i — dz = i Rehan zy ~~~ (12) 
0 Sbacking backing 25packing 


A scaling factor should, therefore, be included in the lateral 
electronic resistance. The modified lateral electronic resistance 
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aS 


ip fr iy 
ii fini 
p 


X, channels and lands, 20 nodes/each 
X, backing layers, 10 nodes/each. 

X, catalyst layers, 3 nodes/each. 

X, membrane, 10 nodes. 

Y, 100 nodes. 

Z, 20 nodes. 


Total : 152,000 nodes. 


Fig. 2. Numerical mesh used for 3D simulations. 


is 
Z 
Renate = Rae 
2ôbacking 
0.5 
=16.7— = -139m2 cm? (13) 


2mm x 0.3mm 


From Eqs. (10) and (13), it is clear that both through-plane and 
in-plane electronic resistances are sufficiently large enough to 
affect the current distribution, not only globally on the cell per- 
formance but also locally between a channel and a land in the 
flow direction. 

The parameter of Zchannel /26backing in Eq. (13), the ratio of the 
half width of the flow channel to the doubled backing layer thick- 
ness, is thus an important variable in the portable DMFC design. 
At a given flow channel width, the in-plane ohmic drop will 
increase when decreasing the backing layer thickness. Numeri- 
cal results presented in the next section reveal that the backing 
layer functions not only to distribute reactants to catalyst layer 
but also to provide a path for lateral electron conduction. Note 
also that the contact resistance causes a resistance increase in 
the through-plane direction and has a minimal effect on the in- 
plane direction. However, the effect of the contact resistance is 
ignored in this study, due to relatively low current density in the 
DMFC. 


3.2. Numerical results and discussions 


In order to solve the model in Fluent, a numerical mesh with 
a total of 152,000 grid points for the 3D single channel DMFC 
domain is generated in Gambit software, as shown in Fig. 2. 
This mesh has been tested to be sufficient by a careful grid- 
independence study. As convergence criteria, we set the residuals 
for all governing equations less than 1076. In addition, selected 
physical properties and electrochemical parameters most rele- 
vant to the present simulation results are listed in Tables 1 and 2, 
respectively. 

As shown in Eq. (13), the in-plane electronic resistance is 
determined by the electronic conductivity, backing layer thick- 
ness and flow channel width. A parametric study is conducted 
to reveal its effect on current density distribution and cell per- 
formance. The variation of the in-plane electronic resistance is 


Table 1 
Three-dimensional cell geometry and operating conditions 
Cell length 0.1m 
Cell width 2x 107° m 
Anode channel width 1x 10m 
Anode backing thickness 300 x 1076 m 
Anode catalyst thickness 10 x 107 m 
Membrane thickness 50 x 10-6 m 
Cathode catalyst thickness 10 x 107 m 
Cathode backing thickness 300 x 10-6 m 
Cathode channel width 1x 107-3m 
Operating temperature 60°C 
Anode channel pressure l atm 
Cathode channel pressure l atm 
Flow stoichiometry of anode channel 2 
Flow stoichiometry of cathode channel 3 
Inlet methanol concentration at anode 2000 mol m~3 (2 M) 
Inlet liquid saturation at anode 100% 
Inlet liquid saturation at cathode 0% (fully humidified air) 
Interfacial liquid saturation at the cathode 0% 

backing layer 
Cell operating voltage 0.4 V 


carried out through changing the backing layer conductivity. The 
effect should be the same through changing the backing layer 
thickness and/or the flow channel width, while the latter must 
involve generation of new numerical meshes, which is rather 
time consuming. Current density distributions in the middle of 
the membrane under three different electronic conductivities of 
the backing layer are shown in Fig. 3. For the electronic con- 
ductivity of 300 S m7}, the maximum current density is located 
under the edges between the gas channel and the current col- 
lecting land, as this location offers the best combination of easy 
access by methanol and a short path for electron transport onto 
the land. It is clearly shown that in-plane electron transport in 
the backing layer plays a critical role in determining the current 
density profile between the channel and the land. 

As the electronic conductivity increases to 1000 S m7!, the 
land effect is weaker, while the methanol transport phenomenon 
in the in-plane direction becomes more dominant. By assum- 
ing an infinitely large electronic conductivity, i.e., completely 
neglecting electron transport, methanol transport becomes the 
sole factor determining the current density distribution; as a 
result, the maximum current density always occurs in the middle 
of the gas channel. The numerical results further indicate that 
the widths of the flow channel and current collecting land are 
key optimization parameters for better cell performance. 

Current density profiles along the in-plane direction in the 
mid-length of the cell under different electronic conductivities 


Table 2 
Electrochemical parameters 


96,485 C mol"! 
8.314 J (mol K)~! 


Faraday constant 
Gas constant 


Anodic transfer coefficient 0.239 

Cathodic transfer coefficient 0.875 

Proton conductivity of membrane 0.1 x 10? S m7! 
Electronic conductivity in current collecting land 20,000 S m7! 
Electronic conductivity in GDL 300-00 S m=! 
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0.170 
0.165 
0.159 
0.154 
0.149 
0.144 
0.138 
0.133 
0.128 
0.123 
0.117 
0.112 
0.107 
0.102 
0.096 
0.091 
0.086 
0.081 
0.075 
0.070 


Opasne = 3008 /m Osecting = 10008 /m 


O backing = °° 


Fig. 3. Current density distribution (A cm~*) in the middle of membrane under 
different electronic conductivities in the backing layer (at 0.4 V). 


of the backing layer are shown in Fig. 4. In order to neglect the 
electron transport effect, the electronic conductivity of the back- 
ing layer has to reach at least 1000 S m™!, at which the relative 
error in the average current density is less than 5% as compared 
to the case with an infinitely large electronic conductivity. The 
average current densities at cell voltage of 0.4 V with different 
electronic conductivities are listed and compared in Table 3. 
The electronic resistance in the backing layer will not only 
affect the current distribution in the width direction, but also that 
in the flow direction. The current density distributions in the mid- 
dle of the membrane under different electronic conductivities in 
the backing layer in the flow direction are presented in Fig. 5. 
The current density distributions are all showing a fast increase 
near the inlet due to the strong methanol crossover, and then 


0.2 


Flow Channel 


0.18 


0.16 


0.14 


0.12 


0.1 


Current Density (A/cm’) 


——- Infinite 


0 0.0005 0.001 


Cell Width (m) 


0.0015 0.002 


Fig. 4. Current density distribution in the width direction in the mid-length of 
the cell for different electronic conductivities in the backing layer. 


Table 3 
Average current densities at cell voltage of 0.4 V with various electronic con- 
ductivities of backing layer 


Electronic Average current Relative error 
conductivity (S m~!) density (A cm~?) (%) 
œo 0.142 - 
1500 0.137 3.5 
1000 0.135 4.9 
500 0.133 6.3 
300 0.132 7.0 


300S/m 
— — = 500S/m 
—-—-— 1000S/m 

1500S/m 
—— — Infinite 


Current Density (A/cm?) 


0 0.025 0.05 
Cell Length (m) 


0.075 0.1 


Fig. 5. Current density distributions in middle of cell width along the flow 
direction for different electronic conductivities in the backing layer. 


Current Density (A/cm?) 


Near Inlet 
— — = Inthe Middle 
—+—-= Near Outlet 


0 0.0005 0.001 


Cell Width (m) 


0.0015 0.002 


Fig. 6. Current density distribution in the middle of the membrane in both flow 
and width directions (backing = 300S m~!). 
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become flat because of weak methanol crossover afterwards. 
Moreover, with lower electronic conductivity, the current density 
in the membrane shifts to a lower level as well. This performance 
reduction is caused by the increasing electrical resistance in the 
backing layer and thus, larger ohmic loss. It is evident that the 
electronic conductivity shows more important effect on the cur- 
rent density distribution than methanol crossover under the given 
conditions. 

Variations in the current distribution in both flow and width 
directions are presented in Fig. 6 for the backing layer electronic 
conductivity of 300S m7! and cell voltage of 0.4 V. In regions 
near the inlet and in the middle of the cell, with higher con- 
ductivity of the bipolar plate than the backing layer, the highest 
local current density occurs at the edges between the land and 
flow channel. On the other hand, near the outlet, due to the severe 
methanol transport limitation in the land area, the highest current 
density shift towards the middle of the flow channel. 


4. Conclusions 


Based on the 3D DMFC model, the electron transport effect 
in the DMFC is numerically explored, with focus on the impact 
of in-plane electron transport on the current distribution and 
cell performance. It was found that the electronic resistance of 
the backing layer, affected by backing layer electronic conduc- 
tivity, backing layer thickness and flow channel width, played 
a relatively important role in determining the current den- 
sity distribution and cell performance. In order to ignore the 
electron transport effect on the average current density, the min- 
imum electronic conductivity of the backing layer has to be 
1000S m7!, with the relative error in the average current den- 


sity less than 5%, under the given conditions. However, such 
effect will become stronger under higher current density condi- 
tions. Future work is needed to examine effects of anisotropic 
electronic conductivity typically found in fuel cell backing layer 
materials. 
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